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Abstract 

Deformation  of  the  wire  used  in  the  solenoidal 
windings  of  an  inertially  confined  pulsed  high  magnetic 
field  generator  is  potentially  the  limiting  factor  for  the 
magnitude  and  duration  of  the  magnetic  field  produced. 
The  rising  magnetic  pressure  at  the  wire  surface  can 
become  large  enough  to  cause  the  cross  section  of  the 
wire  to  deform  on  a  time  scale  shorter  than  overall 
solenoid  disassembly  time.  This  may  result  in  short 
circuiting  due  to  insulator  breakage  and/or  physical 
contact  of  adjacent  windings.  An  analytic  approximation 
modeling  the  deformation  dynamics  is  presented  which 
takes  into  account  both  inertial  and  material  yield  strength 
effects.  The  model  is  validated  by  comparison  to  two 
dimensional  magnetohydrodynamic  simulations  of  the 
process  by  Numerex's  MS  Windows  version  of  AFRL's 
MACH2.  Cases  where  yield  strength  has  a  negligible 
effect  on  the  deformation,  and  where  yield  strength  is 
significant  are  considered. 


wire  is  calculated  assuming  conservation  of  global  energy 
and  momentum. 

The  relationship  between  P  used  in  the  model  and  the 
central  solenoid  magnetic  field  Bl  must  take  into  account 
the  fact  that  the  local  magnetic  field  at  the  wire  surface  is 
higher  than  Bx.  To  determine  an  effective  P,  we  equate  the 
initial  pdV  work  performed  per  winding  of  an  ideal 
solenoid  by  expanding  its  radius  by  dR.  to  the  work 
performed  by  pressure  P  to  our  square  wire  over  the  same 
displacement, 

o 2 

—L  x  2n  RoSdR  =  Px  2nR0AodR  (1) 
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Here,  R0  is  the  initial  solenoid  radius,  S  is  the  wire 
spacing,  /-0  is  the  initial  wire  radius,  A0  is  the  initial  side 
length  of  a  square  wire  of  the  same  area,  and  pQ  =  An  x 
10'7  H  m1.  From  this, 


P  =  A5[ 

Aq  2p0 


where 


A0  =  yFrr0  (2) 


I.  INTRODUCTION 

Consider  a  long  solenoid  wound  from  round  wire  made 
of  metal  with  a  given  yield  strength  and  having  a  diameter 
small  compared  to  the  solenoid  diameter.  We  wish  to 
analytically  model  the  deformation  magnitude  of  the  wire 
due  to  an  internal  magnetic  pressure  high  enough  to  cause 
rapid  disassembly.  We  represent  the  wire  in  planar 
geometry  by  a  square  one  of  the  same  area  that,  by 
assumption,  remains  rectangular  during  deformation  to 
make  the  problem  more  tractable.  The  solenoid  radius  is 
assumed  to  be  large  enough  to  permit  hoop  stress  effects 
to  be  neglected  during  the  time  of  interest.  The  magnetic 
pressure  is  represented  by  a  uniform  pressure  applied  to 
one  side  of  the  wire  that  rises  instantaneously  to  a 
constant  value  P  at  time  t  =  0.  A  free  standing  wire  and 
one  bounded  by  a  rigid  wall  on  the  side  opposite  the 
pressurized  side  are  considered.  These  are  intended  to 
represent  a  free  standing  solenoid  and  one  reinforced  by  a 
surrounding  tube,  respectively.  The  dynamic  width  of  the 


The  models  are  checked  against  the  code  MACH2[1][2] 
for  a  range  of  conditions  to  establish  their  suitability  for 
providing  general  design  criteria  for  similar  geometries. 


II.  DEFORMATION  OF  A  FREE 
RECTANGULAR  WIRE  DUE  TO 
PRESSURE  ON  ONE  SIDE 

Referring  to  “top”  and  “bottom”  as  illustrated  in  Fig.  1, 
we  first  calculate  the  deformation  of  a  free  rectangular 
wire  (for  generality)  due  a  uniform  and  constant  pressure 
P  on  the  top  surface,  assuming  the  wire  remains 
rectangular.  We  define  our  (x,y)  coordinates  to  be  in  the 
accelerating  reference  frame  of  the  wire,  with  the  origin  at 
the  center  of  mass.  The  dynamic  wire  width  is  A  (A0 
initially),  the  height  is  B  (B0  initially),  and  the  center  of 
mass  speed  is  V.  Uniaxial  yield  strength  Y  and  density  p 
are  assumed  constant.  We  further  assume  the  simplest 
possible  expression  for  the  material  velocity  field  in  the 
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center  of  mass  frame  v  consistent  with  our  rectangular 
cross  section  assumption, 

dA  x  dB  y 

"*=dF3  v’  =  StB  '~  =  °  (3) 

From  incompressibility. 
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where  the  integrals  are  over  the  wire  cross  section. 
Conservation  of  energy  implies, 


dW^dEp  dK_ 
dt  dt  dt 
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Substituting  the  v,  D,  S  components  from  Eq.  (3),  Eq.  (6), 
and  Eq.  (5),  respectively,  into  the  energy  terms,  and 
performing  the  integrations,  this  expression  becomes, 


1  „  .  dB  /2"  //I  dA\2  (  1  dB\2 

2  dt  V  3  \J  \  A  dt  J  \B  dt  J 

,  „  rrdV  p  d  \^(dA\2  ( dB \ 2" 
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Figure  1.  Geometry  for  analytic  model  of  free  standing 
solenoid  wire  deformation. 


The  Levy-Mises  equation  for  rigid  plastic  flow[3]  is 
assumed. 


S  =  y  3  YD/ VD  •  D, 


where  S  is  the  deviatoric  stress  tensor  (ST  =  0),  and  D  is 
the  deviatoric  strain  rate  tensor, 


D  =  i(Vv  +  Vvx)-i(V-v)I 


Conservation  of  momentum  implies, 

PA  =  pA0B0~  (14) 

This  allows  us  to  eliminate  the  terms  with  V  in  Eq.  (13).  B 
may  be  eliminated  via  Eq.  (4),  and  the  result  integrated  by 
time  to  obtain. 


ds  j24(P/2-2Y/VS)s  ,  ,  A 

—  =  \  — 77-0 — o - .0  0  v  where  e  =  In  — —  (15) 

dt  V  p  ( />,•;  (.  -  .1,*/  -  )  A0  1 


e  can  be  recognized  as  plastic  strain.  Solving  for  time 
gives, 


Incompressibility  implies, 


V  •  v  =  0 


With  this,  substituting  Eq.  (3)  into  Eq.  (6)  gives, 


p  (Bq  +  Ag ) 


I  Bq6  2s  +  Ao62s 


Y  24  (P/2  -  2Y/>/3)  Jo  V  s{Bl  +  Al) 

(16) 

The  integrand  may  be  approximated  by  the  truncated 
Laurent  expansion  in  s  , 


Dyv  ~  B  dt 


for  the  nonzero  D  terms. 

The  time  derivatives  of  the  external  work,  plastic 
work[3],  and  kinetic  energy  per  unit  length  of  wire  are, 
respectively. 


_i  i  /A  a2\  a  fa  q3\  s 

s  °  '  •  ('-•t)"  ■  (t  t)  ";  (17) 


This  permits  integration  giving  the  following 
approximation  accurate  for  small  strains  (£  -  0. 1 ), 


V  3  P  -  l\  3  V'  [  3  5  V  2 
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where  a  = 


Bl-Aj 
Bl+A  5 


7  V  3  2  J 


where  ^ifr|  ^ 


We  see  from  the  denominator  of  the  root  that  no  flow 

occurs  if  Y>  3  2  P/4.  We  see  from  the  denominator  of  the  root  that  the  yield 

A  =  B0  for  an  initially  square  wire,  and  Eq.  (16)  may  be  strength  must  be  twice  as  high  in  the  supported  case  ( Y  > 
written,  3Vl  P/2)  vs.  the  unsupported  one  to  avoid  flow.  /J  =  3/5  for 


3 P  -  iy/3 Y  Jo 


/cosh  (2  s2)ds  (19) 


And,  with  a  =  0,  Eq.  (18)  becomes, 


3 P  -  1\  3) 


+  ^j  for  .In  B0  (20) 


III.  DEFORMATION  OF  A  BOUNDED 
RECTANGUEAR  WIRE  DUE  TO 
PRESSURE  ON  ONE  SIDE 


an  initially  square  wire. 

IV.  COMPARISON  WITH  MACH2 
SIMULATIONS 
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Figure  2.  Geometry  for  analytic  model  of  externally 

bounded  solenoid  wire  deformation.  Figure  3.  MACH2  results  at  10%  spread  for  two  cases  of 

solenoid  wire  deformation. 

A  similar  analysis  is  performed,  but  for  a  wire  bounded 

from  below  as  shown  in  Fig.  2.  Qualitative  differences  in  MACH2  simulations  run  in  planar  geometry  representing 

the  calculation  are  that  energy  terms  are  evaluated  in  the  a  soienoid  wound  from  wire  of  radius  r0=  3.18  X  10  3  m 

stationary  frame  and  conservation  of  momentum  is  not  and  having  a  gap  between  turns  equal  to  10%  of  the  wire 

needed  to  simplify  the  energy  balance  equation.  The  diameter  are  compared  to  our  analytic  free  and  bounded 

corresponding  time  to  reach  strain  8  in  this  case  is,  wire  models.  We  consider  four  cases  with  B,  =  35  T  and 

Bl  =  50  T  using  half-hard  Cu  and  90%  cold  worked 
Glidcop™  ATO3  dispersion  strengthened  Cu  alloy  AL- 
ds  (21)  15[4][5].  The  latter  has  a  density,  compressive  equation 

of  state,  conductivity,  and  elastic  modulus  similar  to  Cu, 
.  .  but  a  higher  yield  strength.  The  SESAME  Cu  equation  of 

w  ic  is  approximate  or  £  .  y,  state  table[6]  and  Desjarlais'  modified[7]  Lee-More  Cu 

electrical  conductivity  model[8]  are  used  for  both 
materials.  The  Steinberg-Cochran-Guinan  strain  rate 


independent  half-hard  Cu  model[9]  is  used  for  Cu's 
elastic -plastic  behavior.  The  same  S-C-G  shear  modulus 
is  used  for  AL-15  but,  lacking  high  strain  rate  data,  AL- 
15's  uniaxial  yield  strength  is  set  to  its  published  value  of 
3.40  x  108  Pa.  Periodic  boundary  conditions  are  invoked, 
corresponding  to  the  wire  being  one  of  an  infinite  array  of 
wires  with  center-to-center  spacing  of  7  X  10’3  m.  The 
upper  boundary  condition  at  y  =  0. 10  m  is  defined  to  have 
a  magnetic  field  of  B  =  B\X,  keeping  the  same  “up-down” 
convention  as  the  analytic  model.  The  lower  boundary 
condition  is  at  y  =  -0.10  cm  with  B  =  0.  The  simulations 
available  to  date  are  for  free  wires  only.  Fig.  3  illustrates 
the  wire  outlines  for  two  cases  at  the  time  the  wires 
expand  by  10%  and,  therefore,  short  to  the  neighboring 
turns. 

We  compare  the  time  t  required  for  the  wire  to  spread 
by  10%  in  the  different  models.  With  the  help  of  Eq.  (2) 
and  Eq.  (15),  parameters  for  the  analytic  models  common 
to  all  cases  considered  are, 

S  =  7.00  x  1C)-3  m  p  =  8930  kg  nr3 

Bo  =  Ao  =  5.64  x  10  3  m  e  =  0.0953 

(23) 

The  equivalent  pressures  for  B\  =  35  T  and  B\  =  50  T  are, 
from  Eq.  (15),  P  =  6.05  x  10s  Pa  and  P  =  1.24  x  109  Pa, 
respectively.  Using  the  high  strain  rate  initial  yield 
strength  from  half-hard  Cu’s  S-C-G  model  and  the 
published  yield  strength  of  AL-15,  we  use  Y  =  1.23  X  108 
Pa  and  Y  =  3.4  x  108  Pa,  respectively,  for  the  analytic 
models.  Table  1  summarizes  the  spread  times  determined 
by  Eq.  (20),  Eq.  (22),  and  MACH2. 


Table  1.  Times  for  wire  to  spread  10% 


Case 

Free 

Bound 

MACH 

Cu,  35  T 

7.50 

6.85 

7.0 

Al-15,  35  T 

Stable 

10.11 

10.1 

Cu,  50  T 

4.36 

4.46 

4.7 

AL-15,  50  T 

6.34 

5.08 

5.7 

Y  plays  only  a  small  role  in  the  dynamics  and  the  three 
models  are  in  reasonable  agreement  in  all  cases  but  the 
second.  One  may  infer  from  the  small  displacement  of  the 
bottom  of  the  wire  for  Cu  at  50  T  illustrated  in  Fig.  3  that 
a  bounding  wall  does  not  have  enough  time  to  affect  the 
wire  width  by  the  time  the  wire  spreads  10%  in  the 
inertially  confined  limit.  Both  analytic  models  reflect  this 
accurately. 

It  comes  as  a  surprise  in  the  second  case,  however,  that 
MACH2  agrees  much  better  with  the  analytic  bounded 
wire  model  than  with  the  free  wire  model  that  was 
intended  to  represent  it.  In  fact,  the  latter  model  predicts 
that  the  wire  does  not  deform  at  all  in  this  case  ( Y  >  3'/2 
P/4).  It  appears  that  in  cases  such  as  this  where  Y  does 
play  a  significant  role  in  impeding  expansion,  Eq.  (22) 
provides  the  best  guidance.  Though  not  yet  fully 
understood,  Fig.  3  suggests  an  explanation.  The  wire's 
shape  at  10%  expansion  varies  little  with  the  relative 


values  of  Y  and  By.  Furthermore,  the  deformation  occurs 
primarily  in  the  top  half  before  the  bottom  is  substantially 
displaced,  at  least  for  the  cases  explored  here.  Placing  a 
lower  boundary  on  the  analytic  model  apparently  offsets 
the  error  introduced  by  its  unrealistically  homogeneous 
deformation.  It  is  not  clear,  however,  that  this  would 
continue  to  be  the  case  for  significantly  greater 
deformations,  or  if  P  were  much  closer  to  the  critical 
value  for  flow  in  the  bounded  model  (3.93  x  108  Pa  for 
AL-15).  These  cases  will  be  addressed  in  future 
calculations  to  better  characterize  the  range  of  validity  of 
Eq.  (22). 


V.  REFERENCES 

[1]  M.  H.  Frese,  “MACFI2:  A  two-dimensional  magneto¬ 
hydrodynamic  simulation  code  for  complex  experimental 
configurations”,  Mission  Research  Corp.,  Albuquerque, 
NM,  Tech.  Rep.  AMRC-R-874,  1987. 

[2]  R.  E.  Peterkin,  Jr.  and  M.  H.  Frese,  “A  material 
strength  capability  for  MACH2”,  Mission  Research  Corp, 
Albuquerque,  NM,  Tech.  Rep.  MRC/ABQ-R-1 191,  Nov., 
1989. 

[3]  A.  S.  Kahn  and  S.  Huang,  Continuum  Theory  of 
Plasticity.  NY:  John  Wiley  &  Sons,  1995. 

[4]  T.  J.  Miller,  S.  J.  Zinkle,  and  B.  A.  Chin,  “Strength 
and  fatigue  of  dispersion  strengthened  copper”,  J.  Nuclear 
Mater.,  179-181,  pp.  263-266,  1991. 

[5]  T.  S.  Srivatsan,  N.  Narendra,  and  J.  D.  Troxell, 
“Tensile  deformation  and  fracture  behavior  of  an  oxide 
dispersion  strengthened  copper  alloy”,  Materials  and 
Design},  vol.  21,  pp.  191-198,  2000. 

[6]  S.  P.  Lyon  and  J.  D.  Johnson  (eds.),  “SESAME:  the 
Los  Alamos  National  Laboratory  equation  of  state 
database”,  LANL,  Los  Alamos,  NM,  Tech.  Rep.  LA-UR- 
92-3407,  1992. 

[7]  S.  E.  Rosenthal,  M.  P.  Desjarlais,  R.  B.  Spielman,  W. 
A.  Stygar,  J.  R.  Asay,  M.  R.  Douglas,  M.  H.  Frese,  R.  L. 
Morse,  and  D.  B.  Reisman,  “MHD  modeling  of 
conductors  at  ultrahigh  current  density”,  IEEE  trans. 
Plasma  Sci.,  vol.28,  pp.  1427-1433,  2000. 

[8]  Y.  T.  Lee  and  R.  M.  More,  “An  electron  conductivity 
model  for  dense  plasma”,  Phys.  Fluids,  vol.  27,  pp.  1273- 
1286,  1984. 

[9]  D.  J.  Steinberg,  S.  G.  Cochran,  and  M.  W.  Guinan, 
“A  constitutive  model  for  metals  applicable  at  high-strain 
rate”,  J.  Appl.  Phys.,  vol  51,  pp.  1498-1504,  1980. 

[10]  D.  J.  Steinberg,  “Equation  of  state  and  strength 
properties  of  selected  materials  (revised)”  LLNL, 
Livermore,  CA,  Tech.  Rep.  UCRL-MA-106439,  change 
1,  Feb  1996. 


